Coulomb blockade of a noisy metallic box: A realization of Bose-Fermi Kondo models 
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We focus on a metallic quantum dot coupled to a reservoir of electrons through a single-mode point 
contact and capacitively connected to a back gate, by including that the gate voltage can exhibit 
noise; this will occur when connecting the gate lead to a transmission line with a finite impedance. 
The voltage fluctuations at the back gate can be described through a Caldeira-Leggett model of 
harmonic oscillators. For weak tunneling between the lead and the dot, exploiting the anisotropic 
Bose-Fermi spin model, we show that zero-point fluctuations of the environment can markedly alter 
the Matveev Kondo fixed point leading to an amplification of the charge quantization phenomenon. 



PACS numbers: 73.23.Hk, 42.50.Lc, 72.15.Qm 
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The single-electron box, the simplest single-electron 
circuit, is probably the ideal system in which to test 
the theory of quantum charge fluctuations. The box has 
been vividly investigated theoretically because this is a 
model system for understanding electron-electron inter- 
actions and because the quantum fluctuations in the box 
are analogous to the Kondo effect [J- In practice, charge 
fluctuation measurements necessitate a (large) dot at the 
micron scale with a very dense spectrum |2j; This is re- 
ferred to as a metallic dot. In earlier theoretical treat- 
ments of the metallic dot coupled to a reservoir leadHHl, 
the gate voltage was treated as a fixed parameter of the 
Hamiltonian. As depicted in Fig. 1, in this Letter we 
consider the more general situation where the source of 
the gate voltage is placed in series with an impedance 
Z(cj) of the gate lead showing resistive behavior at low 
frequency. From the fluctuation-dissipation theorem this 
resistance will introduce noise in the gate voltage even 
at zero temperature. Below, dot and lead are coupled 
through a single- mode quantum point contact (QPC). 

Since the charge on a quantum box can now be mea- 
sured up to few thousandths of an electron^, it is 
then important to ask to what extent zero-point fluctu- 
ations of the electrical environment affects the (Kondo) 
ground state energy, i.e., the capacitance, of a metallic 
dot already coupled to a reservoir of electrons through 
a single-mode QPCyLj. While the transport through a 
single ultrasmall tunnnel junction coupled to an electro- 
magnetic (dissipative) environment has been thoroughly 
investigated 5], to our knowledge a discussion for the 
Coulomb blockade oscillations of a metallic box subject 
to a fluctuating gate voltage seems not to be available. 
According to conventional wisdom, tunneling rates of 
quasiparticles become strongly suppressed in the presence 
of dissipative environments because an additional energy 
is needed to excite the environment]^ ; Hence, this should 
certainly modify the Kondo picture ascribed by Matveev 
for electron tunneling between the island and the leadQ]. 

The back gate here is embodied by a voltage fluctuat- 
ing around its average value V g . We treat the full dynam- 
ics of the fluctuations in a self-consistent manner; SV g (t) 
representing the voltage fluctuations (or the charge dis- 
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FIG. 1: Circuit diagram of the metallic dot coupled to a 
bulk lead and subject to a fluctuating back-gate voltage 
Vg + SV g (t). The circuit is supposed to exhibit an Ohmic resis- 
tance R = Z(u> = 0) . We incorporate the finite resistance in a 
Hamiltonian fashion, through a long dissipative transmission 
line with capacitance Ct and inductance L: R = (L/Ct) 1 ^ 2 ■ 

-placement up to a factor) on the back gate is assumed 
to be given by the sum over the normal coordinates of a 
Caldeira-Leggett type bath of harmonic oscillators 7J : 
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This mapping is entirely justified in the context 
of a long dissipative transmission line; The precise 
correspondance^SJ between mechanical and electrical 
quantities can be found, e.g., in the recent Ref. H In 
the case of a linear circuit element with Ohmic resis- 
tance R, one expects Johnson-Nyquist correlation func- 
tions (SVgSVg) u fa URjj coth {huj/2k B T)\ This usually 
refers to as an Ohmic dissipative environment. For con- 
venience, we have extracted the dimensionless parame- 
ter 7 = R/ Rk which naturally characterizes the cou- 
pling between the bosonic bath and electrons of the 
box; Rk = h/e 2 ps 25.8fcS7 being the quantum of re- 
sistance. Having in mind a long dissipative transmission 
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line with capacitance Ct and inductance L, the Ohmic 
resistance is R = y/L/Ct- Furthermore, the couplings 
Xj, the mass and the frequencies of the oscillators en- 
ter only through the bath's spectral density J(w) = 
7 2 - Xj 2 TrS(uj — u)j)/(2Mcjj) = TiRlu in the quantum 

limit; cu c = y/l/(LCt) is a natural low-frequency cutoff 
for the transmission line. Using the circuit of Fig. 1, the 
Coulomb term for the quantum box takes the form: 

H c = E C (Q - N) 2 - E C N 2 - eQSV g . (3) 

Here, eN — V g C g is proportional to the mean gate volt- 
age and eQ depicts the charge on the dot. Here, the 
charging energy of the granule takes the form E c — 
e 2 /(2Cx;) where Cs = C g + Cf, C g denotes the classi- 
cal capacitance between the back-gate and the grain and 
Ci is the capacitance between the lead and the island. 

Let us first study the case of a low-transparency tunnel 
barrier and focus on the point N=l/2 where the states 
with Q=0 and Q=l are energy degenerate, and charge 
fluctuations in principle can be large. We introduce the 
small parameter h — eV g — e 2 /(2C g ) oc (N — 1/2) to 
measure deviations of N from the degeneracy point. Our 
treatment of electron tunneling is based on the assump- 
tion that the tunneling Hamiltonian takes the formal 



H t = J2\ t \c t dk ci q S + + h.c, 

k,q 



(4) 



where |£|c^ fe Ci g is the tunneling term transferring an elec- 
tron from the lead (1) to the (large) dot (d), and where S + 
is an operator raising the charge Q on the dot from to 1. 
The tunneling matrix element t is assumed not to depend 
on the momenta k and q. Below, we will discuss in detail 
the case of spinless electrons, assuming that a strong in- 
plane magnetic field has been applied, and as mentioned 
previously the case of a point contact most probably with 
one conducting transverse mode[l|, OH. We are led to the 
obvious identification Q = S z + 1/2, and finally to 

H c = -hS z - yS z $; (5) 

^(*) — e J2j XjXj(t) now stands for the bosonic vari- 
able. It is also appropriate to visualize H t as a transverse 
Kondo Hamiltonian H t = ( J±/2) (s + S~ + h.c.) where 
J_L = 2\t\ and s + — J^k q c l q c dk- Bear in mind that here 
a Kondo flip process means physically the transfer of an 
electron from lead to dot or vice- versa, the effective spin 
index a = I, d is in fact the position of an electron in the 

structure, and then s z = J2k, q ( c L c '<z - c dk c dg) ■ 

Hence, in the limit of weak tunneling (\t\p) <C 1, our 
system is embodied by the Bose-Fermi spin Hamiltonian 



H = E H Kin ~ hS Z + H B 
\a=l,d J 

+ ^ (s + S- + h.c.) - 7 $S Z ; 



(6) 



Hft in stand for the usual kinetic terms in the lead and 
dot. The densities of states in the dot and lead have been 
assumed to be equal for simplicity, and are denoted p. It 
is suitable to recall that the Bose-Fermi Kondo model has 
been analyzed in detail in the context of quan tum critical 
points of certain heavy fermion materials [ill llli > but 
to our knowledge Eq. © would constitute its first real- 
ization in mesoscopic structures. We get the following 
Renormalization Group (RG) equations akin to those of 
Zhu and Si [13 and Zarand and Dernier (for e = 0) 
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I = ln(A Q /A) denotes the scaling variable with A D w 
min{_E c , Hlu c } the initial value of the energy cutoff and 
A Q = E c when R = [L = 0). Our bare values obey 
X± = J± P = 2\t\p, X z = J zP = 0, 7 2 = R/R K < 1; 
J z s z S z refers to the induced Ising part of the Fermi 
Kondo coupling^ , and v which is of the order of 1 refers 
to the normalization constant occurring in the boson time 
propagator. Note that here there is no generated trans- 
verse Kondo part for the bosonic coupling. It is obvi- 
ous that since the bosonic heat bath adjusts itself to a 
given (pseudo-)spin configuration this tends to suppress 
the renormalization of the Kondo process, i.e., the coher- 
ent tunneling of electrons from lead to dot. The compe- 
tition between the bosonic and fermionic fields then will 
give rise to two stable quantum phases corresponding to a 
Kondo Fermi-liquid phase (where the Kondo energy scale 
will depend on R) and to a merely bosonic Ising regime. 

It is advantageous to rewrite the RG equations above 
in terms of the two variables and X z — (X z — (v/2)j 2 ). 
As long as 7 2 -C 1, this enables us to recover the 
renowned Kosterlitz-Thouless equation flow: dX±/dl = 
X±X Z and dX z /dl = X± 2 . Hence, this strongly validates 
the idea of a critical value for the external resistance, 
namely i? c [|t|p] w 4|t|pi?x/^, at which the Kondo physics 
will be washed out. Let us underline that this formula is 
only strictly applicable in the case of a low-transparency 
tunnel barrier where iZ c [|t|p] <C Rk, i.e., in the scope 
of validity of the perturbative RG. For R < 4\t\pRK /v, 
since A_l and X z will both grow under renormalization 
a Kondo phase with restoration of the SU(2) symme- 
try emerges. The correction to the classical capacitance 
C = d(Q)/dV g is proportional to the impurity suscepti- 
bility x = d(S z ) I dh. An important point is that when 
applying a strong magnetic field the fixed point cor- 
responds to the usual one-channel Kondo Fermi liquid 
where \ — (TkIR])^ 1 , and the Kondo energy scale can 
be approximated by the one of the completely symmet- 
rical model T K {R] ~ A exp[-l/(A_L - ^7 2 /2)] < E c . 
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FIG. 2: Average dot's charge (Q) versus the rescaled voltage 
N — VgCg/e for weak tunneling (\t\p) <C 1 and in the case 
of a single-mode spin-polarized point contact. When R will 
exceed 7? c [|f|p] « 4\t\pR,K /v, we predict the restoration of 
perfect Coulomb steps|l4j|. (inset) Evolution of the Coulomb 
blockade oscillations for R — and various dot-lead couplings. 

Close to the degeneracy point N = 1/2, the aver- 
age dot's charge (Q) then varies linearly with the ap- 
plied gate voltage V g at zero temperature and the slope 
C = d(Q)/dV g will progressively increase with R (see 
Fig. 2). When R > 4\t\pRx/v, the Kondo temperature 
vanishes which means that \±(l) — > and that 7 re- 
mains almost unrenormalized. Since ($) = 0, we expect 
a perfect jump in (Q) or (S z ) similar to a free spin in 
magnetic field: (Q) = 1 if h > (N > 1/2) and (Q) = 
if h < (N < 1/2). This shows that a reasonably large 
resistance in the electrical circuit can markedly restore 
a perfect Coulomb staircase behavior by suppressing the 
virtual tunneling of electrons between the lead and the 
metallic dot. Note that an analogous phenomenon has 
been reported in the different context of a small "noisy" 
dot (embodied by the lowest unoccupied electron level as 
opposed to our many-body grain) embedded in a small 
mesoscopic ring (described by the topmost occupied elec- 
tron level as opposed to our semi-infinite lead) |9J . In the 
case of spinful fermions, i.e., without magnetic field, the 
RG flow would be quite similar to that of Eq. J7J • How- 
ever, the effect of the external impedance should be much 
less spectacular because already in the Kondo phase there 
is a blatant divergence in the capacitance for N ss 1/2 lj 

emerging clearly from the two-channel Kondo model (the 
two channels being the spin polarizations of an electron) . 

A maybe more physical way of understanding the effect 
of the electrical environment is to absorb 5V n in the tun- 
neling term using a unitary transformation 15] following 
Refs. |5|, |6|. In our situation, H t naturally turns into: 

H t = J ± e^8 + S- +h.c, (9) 



where the phase ip(t) — (e/h) f_ SV g (t')dt' at the time 
t represents the integral over the voltage fluctuations felt 
by the dotQ. It becomes (more) transparent that a trans- 
ferred charge onto the dot must be accompanied by zero- 
point fluctuations ( excitations ) of the environment. In 
this formulation of the tunneling Hamiltonian, the Kondo 
coupling J z will be renormalized by the contribution 

J ± 2 J dr sgn{ T ) G(t) e K{T \ (10) 

G(t) = 1/t corresponds to the usual fermionic imaginary 
time propagator. The main difference with the original 
Kondo problem is the occurrence of the phase correlator 
K(t) — (<p(T)(p(0) — p(0) 2 ). An exact calculation of K (r) 
can be performed at zero temperature for the transmis- 
sion line, i.e., for the bath of harmonic oscillators, and for 
large r = it we find K(t) = -2{R/R K ) In (w c |t|)@ ®. 
When R is negligible, K(t) « 0, one easily recovers 
the second line of Eq. JJJ) and then the expected one- 
channel Kondo fixed point. By increasing R, the integral 
Ih/A d T r _2a_1 involved in the above equation (|10|l will 
remain very small whatever the energy A. Hence, this 
leads to the same conclusion as previously: J z (l) (J±(l)) 
will persist to be insignificant up to zero temperature and 
then the (pseudo-)spin S z will satisfy (S z ) = ±1/2. This 
can be entirely attributed to the fact that, for large R, 
the probablity P(t) = e K( - r " 1 to observe zero-point fluc- 
tuations in the environment rapidly decreases for a long 
time t^EI and as a result long-time (low-energy) Kondo 
(tunneling) processes will be eliminated by the noise. 

Now, we shall concentrate on the opposite limit, i.e., 
close to perfect transmission. Here, electronic wave pack- 
ets are clearly spread out between the bulk lead and the 
granule, and hence one could anticipate that the zero- 
point fluctuations of the environment will have more 
difficulty to affect the charge fluctuations of the quan- 
tum box. To make it more quantitative, we can resort 
to bosonization and tackle the one-dimensional problem 
emerging at the QPC along the lines of Refs. Ill llOL Con- 
sequently, the Coulomb term becomes 

H c = ^ (0(0) - V^N) 2 - E C N 2 - A^WWv (11) 

where the boson e</)(0)/v / 7r describes the charge on the 
granule and the coordinate x = refers to the entrance 
of the dot (or QPC). First, we already note that exactly 
at perfect transmission (p\t\ = 1), as long as (SV g ) = 0, 
the capacitive coupling with the back gate will definitely 
impose (0(0)) w t/ttN, and then the grain's charge will 
vary linearly with the (mean) back-gate voltage whatever 
the external resistance (inset in Fig. 2). Obviously, the 
previous qubit basis S z = ±1/2 would not be appropriate 
in that limit. The weak backscattering contribution at 
the QPC then can be manipulated as|l| 

H bs =- — \r\ cos f V^56(0)) cos (2?rJV) , (12) 
wa V / 
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a is a short-distance cutoff, v F denotes the Fermi ve- 
locity, the small reflection amplitude obeys |r| = (1 — 
(pltl) 2 ) 1 / 2 <C 1, and the charge fluctuation field is 
6(f)(0) = 0(0) — t/ttN. For spinless fermions, the (first- 
order) correction to the ground state energy naturally 
takes the form SE 1 = -v F / '(7ra)Me~ M,50(O)2> cos {2nN). 

Let us first recall that when R/Rk — > 0, (6(f)(0) 2 ) 
is large but finite, and a simple calculation provides 
(6(f)(0) 2 ) = -1/(2tt) In (apE c /v F ) where \i = e c and C « 
0.5772 is the Euler's constant; the phonons with energies 
below E c are pinned down by the interaction term. The 
periodic correction to the ground-state energy then reads 

M 

8 Ei = -—hE c cos{2ttN). (13) 

The average charge in the dot is given by e(Q) — eN — 
e/(2E c )d6E 1 /dN = eN-fie\r\ sin (2?riV). One recovers a 
curve similar to the dashed line in Fig. 2 even though the 
amplitude of the charge oscillations is much smaller (i.e., 
this is much closer to the dotted straight line). Now, we 
revisit the quantum fluctuations of the charge displace- 
ment for larger external impedances. The part of the 
Coulomb term containing 6(f)(0) can be turned into 

H c = ^6(f>(0) 2 ~ e<5<K0) Y>Xj, (14) 

3 

with f3j = 'yXj / s/tt. It is suitable to note the resemblance 
with a Brownian particle in a harmonic potential. By ex- 
ploiting this analogy0,0|, we can check that the follow- 
ing correlation function (^(0,4)5^(0,0)) cx -R/ (E c 2 t 2 ) 
almost vanishes for long times and hence that the noise 
is rather irrelevant for |r| <C 1. The charge displace- 
ment (6(f)(0) 2 ) will be still dominated by the logarithmic 
phonon contribution above cutoff at the time h/E c . Zero- 
point fluctuations of the transmission line cannot funda- 
mentally affect the result of Eq. (|13f) . but for very large 
and probably unrealistic impedances R/Rk —> +oo, 
we could expect that the logarithmic contribution of 
(6(f)(0) 2 ) will be cutoff at r c < %/E c . In the absence 
of the (Coulomb) term E c 6(f>(0) /n, for long times one 
would still obtain a logarithmic growth of the 6(f)(0) cor- 
relation function^, (S(f>(0,t)S(j)(0,0)) cx -(R K /R) Int. 

More physically, the scattering of electrons at the QPC 
produces a phase shift in the wavefunctions ipk{ x ) = 
cos(fc|x| —6) which is naturally related to the average 
charge entering onto the dot through the Friedel sum 
rule 6 — ir(Q)/e. As long as \r\ <C 1, this means that 
we can approximate (Q) / e ~ N whatever the exter- 
nal Ohmic resistance (again (6N) = 0). This ensures 
the following Friedel oscillation of the electron density 
p(x) ps E c /(hv F ) cos(2k F \x\ — 26) for |x| <C v F /E c with 
26 ~ 2irN; the restriction of wave vectors k to the range 
E c /(hv F ) around the Fermi momentum k F is because 
only for those wavevectors elastic reflection from the con- 
tact takes place. We can then infer that the shift in the 



ground state energy — J dx p(x)V(x)8(x), where the po- 
tential V — \r\hv F , (roughly) is still equal to 5E\. 

In closing, we have shown that in the context of a 
metallic granule and a low-transparency tunnel barrier, 
zero-point fluctuations of a transmission line can alter 
the quantum tunneling of electrons at low energy and 
then restore a perfect Coulomb staircase behavior. A 
clear observation of such a suppression of the quantum 
charge fluctuations would necessitate the application of 
a strong in-plane magnetic field, a reasonable external 
impedance (a non- negligible fraction of Rk), and again 
a small transmission coefficient between the metal lead 
and the grain. When approaching the perfect transmis- 
sion, quantum charge fluctuations become so prominent 
such that the electrical environment is rather inefficient. 
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